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FINAL REPORT ON DOPPLER LIDAR RESEARCH 

CONTRACT NAS8-37585 

April 1, 1989 -- September 30, 1992 

DEVELOPMENT OF A GLOBAL BACKSCATTER MODEL 
FOR NASA’S LASER ATMOSPHERIC WIND SOUNDER 


I. RESEARCH OVERVIEW 

A. Summary 

During the Contract Period April 1, 1989, to September 30, 1992, the Earth Systems 
Science Laboratory (ESSL) in the Research Institute at the University of Alabama in Huntsville 
(UAH) conducted a program of basic research on atmospheric backscatter characteristics, leading 
to the development of a global backscatter model. The ESSL research effort was carried out in 
conjunction with the Earth System Observing Branch (ES43) at the National Aeronautics and 
Space Administration (NASA) Marshall Space Flight Center, as part of NASA Contract NAS8- 
37585 under the Atmospheric Dynamics Program at NASA Headquarters. This research 
provided important inputs to NASA's GLObal Backscatter Experiment (GLOBE) program, 
especially in the understanding of global aerosol life cycles, and to NASA’s Doppler Lidar 
research program, especially the development program for their prospective space-based Laser 
Atmospheric Wind Sounder (LAWS). 

B. Personnel 

Listed below are the names, positions, and tenures of the ESSL personnel who were 
involved in this research effort. 

Scientists 

David A. Bowdle 

Maurice A. Jarzembski 

Dean R. Cutten 

Database Managers 

Steven F. Williams 

Michael Moore 

Helen Conover 


Research Scientist 
Scientist 

Research Scientist 


April 1989 
April 1989 
June 1990 


September 1992 
December 1989 
September 1992 


Research Associate 
Student Assistant 
Research Associate 


April 1989 
July 1990 
March 1991 


March 1991 
March 1991 
September 1991 


Staff As sistants (Meetings Coordination. Publications. Administrative) 


Dbde Upton 
Laurie S. Collins 
Reena Pearson 
Cindy Taylor 


April 1989 - June 1989 

July 1989 — September 1992 

November 1990— December 1990 
April 1990 — September 1992 



C. Objective and Scope 

The primary objective of this research is to develop a model of the volume backscatter 
coefficient m'^ sr^) for trace atmospheric particulate matter over a wide range of spatial, 
temporal, and spectral scales. The model deals primarily with atmospheric aerosols and clouds, 
with lesser attention to gaseous molecular backscatter. It covers the troposphere and lower 
stratosphere, with special emphasis on low P values for so-called background aerosol populations 
in the middle and upper troposphere. Model-related studies deal with particulate distributions 
from the turbulent sub-microscale (~1 mm to ~10 m) to the global scale (>10,000 km), although 
actual model products span a much more limited range from the synoptic scale (-1000 km) to the 
global scale. Similarly, model inputs cover time periods from the sub-microscale (-1 ps) to the 
climatic (1-10 yr), while model outputs cover the synoptic (-1-10 da) to the climatic temporal 
scales. Research was conducted on aerosol backscatter, extinction, and related optical properties 
from -0.3-15 pm wavelength, but again model products concentrated on Rvalues at selected laser 
wavelengths in the near and middle infrared. 

Development of this model required extensive compilation, careftil review, and creative 
synthesis of aerosol data from a wide variety of sensors and sensor platforms. The deliverable 
products from this research include the backscatter model data inputs, the backscatter model data 
products, the analytical methods used to synthesize the model products from the model inputs, 
and the scientific rationale used to design the analytical methods. 

The ESSL model development effort is an integral part of NASA's GLOBE program. 
GLOBE focused on two major airborne field programs, the Fall 1989 (GLOBE I) and Spring 
1990 (GLOBE II) backscatter survey missions over the Pacific Ocean on a DC-8 aircraft from 
NASA’s Ames Research Center (ARC). These missions included two pulsed lidars, to probe 
aerosol backscatter above and below the aircraft; two continuous-wave (CW) focused Doppler 
lidars, to characterize backscatter at flight altitude; and several auxiliary sensors, to measure 
other physical, chemical, and optical properties of the aerosols at flight level. The instrument 
complement and the theater of operations were chosen to improve NASA's understanding of the 
physicochemical properties and the spatial/temporal distribution of tropospheric aerosols over the 
remote oceans. 

D. Applications 

Backscatter data products from the ESSL model are making important contributions to 
several facets of the LAWS program. In particular, the LAWS project uses ESSL backscatter 
predictions to obtain parametric estimates of LAWS performance for programmatic and 
engineering trade studies on LAWS. ESSL backscatter model products are also used in detailed 
Observing System Simulation Experiments (OSSE's) for LAWS, in order to assess the scientific 
impact of the decisions from these trade studies. 

The results of this research are applicable to similar design studies and performance 
simulations for any earth-based or space-based lidar system that uses atmospheric aerosols as 
passive scattering targets for the measurement of winds or other primary atmospheric variables. 
This research is also providing important insights about the physical, chemical, and optical 
properties, the spatial/temporal structure, and the overall dynamics of the regional- to global-scale 
aerosol system. 



n. RESEARCH DESCRIPTION 
A. Tasks 

ESSL participation in the GLOBE research program covered six general areas. 

1. GLOBE/LAWS Program Leadership ESSL provided expert scientific leadership for 
NASA's GLOBE program and expert scientific representation to NASA's LAWS development 
program. Under this task, ESSL continually provided NASA MSEC and NASA Headquarters 
with state-of-the-art scientific information on atmospheric aerosol systems; their physical, 
chemical, and optical properties; their spatial and temporal distribution; and their probable 
effects on LAWS performance. This information was used in conjunction with NASA MSEC and 
NASA Headquarters to develop a scientifically sound rationale for the overall GLOBE research 
effort and the various elements of GLOBE; to develop realistic research programs based on that 
rationale; and to evaluate or modify those programs based on the evolving needs of the LAWS 
program. ESSL supplemented these efforts by coordinating, admim'stering, and documenting 
periodic meetings of the GLOBE and LAWS Science Teams. 

2. GLOBE Science Team Coordination ESSL coordinated the research efforts of the 
GLOBE Science Team. This activity included designing and directing major GLOBE field 
programs; directing GLOBE data processing and data analysis; coordinating GLOBE special 
sessions at open scientific conferences; and serving as Lead Guest Editor for a special GLOBE 
section of a prominent earth science research journal. 

3. GLOBE Database Development ESSL designed, developed, implemented, and 
maintained a centralized, remotely accessible, interactive GLOBE database at MSEC. This effort 
covered the hardware, data archival software, and data visualization software, as well as data 
compilation, cataloging, archival, and distribution. 

4. GLOBE Model Development ESSL used the comprehensive GLOBE database to 
develop an empirical global-scale model of aerosol backscatter for input to the LAWS 
development effort. This activity involved developing customized data visualization software 
modules; reviewing and quality controlling GLOBE data; intercomparing contemporaneous 
GLOBE data sets; characterizing the physical, chemical, and optical properties of atmospheric 
aerosols firom the GLOBE data; and synthesizing the combined results of the GLOBE program 
into a global-scale aerosol backscatter model. 

5. Aerosol Laboratory Development ESSL developed a new aerosol analysis laboratory 
at MSEC. This activity included designing and implementing the laboratory layout, installing 
MSEC lidars in the laboratory, equipping the laboratory with aerosol analysis accessories, using 
the lidars to measure aerosol optical properties, and developing computer codes that simulate the 
lidar measurement process. 

6. Publications. Presentations. Reviews ESSL presented numerous papers on GLOBE 
results, published several GLOBE papers, and reviewed a variety of proposals and journal 
manuscripts. 

B. Key Results 

1. GLOBE/LAWS PROGRAM LEADERSHIP 

a. GLOBE Program Design D. Bowdle developed the scientific rationale for the entire 
GLOBE research program, and particularly for the GLOBE survey missions. This rationale was 
based on the premise that GLOBE should concentrate on aerosol targets near the detection 



threshold for the current or prospective LAWS design. Using this premise, GLOBE studies for 
the initial LAWS instrument were designed to provide information on aerosol backscatter 
properties at 9.1 pm in the middle and upper troposphere, especially over the remote oceans. 
During late 1991 and all of 1992, the focus of the LAWS program shifted toward smaller, less 
complex, less risky, and less costly versions of LAWS, possibly at a shorter wavelength, such as 
2.1 pm. Using the same premise, GLOBE studies for LAWS began to emphasize higher 
backscatter targets, such as thin clouds, continental aerosol plumes, convective boundary layer 
plumes, and the continental and marine boundary layers, at a variety of wavelengths. This change 
in emphasis could have a profound effect on the design of future GLOBE measurement programs. 
During this same time, D. Bowdle and D. Cutten made major contributions to a new MSFC 
proposal for continuing GLOBE aerosol research by the ES43 Aerosol/Lidar Group. 

*>• GLOBE Science Team Meeting s The GLOBE Science Working Group (SWG) met in 
May 1989, to evaluate readiness for the Fall 1989 GLOBE Survey mission, and in March 1990, to 
review the results of the Fall Survey mission and evaluate readiness for the Spring 1990 Survey 
mission. The SWG met again in May 1991 and April 1992 to review mature results from the 
combined missions. The 1991 meeting featured high-end desktop computers for on-site 
intercomparisons of the measurements from various GLOBE instruments. ESSL organized the 
technical sessions and handled the logistical arrangements for all of these meetings, then published 
the proceedings for most of the meetings. One or more ESSL scientists attended these meetings 
and presented the results of their GLOBE research, emphasiang the typical values and the spatial 
and temporal variability of atmospheric aerosol backscatter. 

c. LAWS Science Team Participation The LAWS Science Team met twice each in 1989, 
1990, 1991, and 1992. ESSL hosted and arranged most of these meetings and published 
Proceedings for all of the meetings. One or more ESSL scientists attended each meeting. 

D. Bowdle joined the Science Team as an Associate Member in early 1989, at the 
invitation of NASA and the LAWS Science Team Leader. Shortly thereafter, he became 
Chairman of the LAWS Backscatter Subcommittee. In these capacities, he provided long-term 
scientific liaison between the GLOBE research program and the LAWS development program. 
At each Science Team meeting, he presented GLOBE status, plans, and recent results, 
emphasiMg implications for LAWS design and performance. He led LAWS Backscatter 
Committee meetings, to determine the current and projected needs of the LAWS program for 
GLOBE data inputs. He also presented Committee findings and recommendations to the Science 
Team, and documented those results for the Meeting Proceedings. 

d. LAWS Design Review D. Bowdle attended numerous quarterly review meetings for 
the LAWS Phase I/n prime contractors. He also attended various presentations on alternate 
LAWS design concepts by other contractors, and prepared written technical reviews for NASA 
on some of these presentations. He also participated extensively in the effort to develop suitable 
scieiitific and engineering design concepts for a descoped version of LAWS. In this capacity, he 
provided information about aerosol backscatter properties and other atmospheric optical 
parameters for a wide variety of laser wavelengths and atmospheric conditions. He also prepared 
a preliminary model of the wavelength dependence and frequency of occurrence for backscatter 
coefiicients and LAWS signal-to-noise ratios at several LAWS power levels and the two prime 
laser wavelengths, 9.1 pm and 2. 1 pm in the dominant tropospheric aerosol types. He presented 
these results, along with an overview of the GLOBE results, to Dr. Shelby Tilford, Director of the 
Earth Science and Applications Division at NASA Headquarters. He also attended the 2 pm 



Solid State Laser Technology Assessment for Remote Sensing of Winds, again at NASA 
Headquarters. 

2. GLOBE SCIENCE TEAM COORDINATION 

a. Field Programs ESSL designed, coordinated, and participated heavily in the GLOBE I 
and GLOBE II airborne backscatter survey missions over the Pacific Ocean. The GLOBE I 
deployment was delayed and re-routed due to instrument problems. Once deployment began, 
most instruments performed well and flight operations were conducted essentially as planned 
under the modified routing. The mission included both long-distance, high-altitude transit flights 
around the Pacific basin, and detailed local flights at selected sites for comparison with 
independent aerosol research stations. A wide variety of meteorological conditions was en- 
countered during the deployment, with corresponding changes evident in aerosol properties. 
After major improvements in all instruments, the Spring survey mission deployed on schedule. 
The general routing, the individual flight plans, and the meteorolo^cal variety for this mission 
were similar to those in the Fall mission, except for a seasonal reversal in hemispheric aerosol 
properties. All instruments performed very well. 

D. Bowdle served as Mission Science Coordinator for both survey missions. Before 
deployment, he developed the scientific framework for each mission; designed flight plan 
concepts to meet key scientific objectives; defined requirements for supporting measurements; 
and assisted in preparing the Mssion operations plan documents. During deployment, he 
collected weather forecast data; participated in preflight planning; provided in-flight 

meteorological observations; directed in-flight sample collection and documentation; contributed 
to in-flight maneuver planning; coordinated near-real-time and post-flight data review; and 
maintained liaison between the Mssion Scientist and the Principal Investigators. After 
deployment, he coordinated data processing, analysis, and validation efforts, and organized 
overall mission debriefings. 

S. Williams served as Mssion Logistics Coordinator and Mission Data Base Manager. 
Before deployment, he developed detailed plans for implementing flight plan concepts; repre- 
sented MSFC on a sdentific, logistical, and political site survey tour, traveled fi^equently to 
NASA ARC and NASA Headquarters for mission planning discussions; monitored interagency 
and international arrangements; arranged for supporting meteorolo^cal and aerosol 

measurements fi'ora numerous ground-based, airborne, and satellite stations; and prepared 
Mission Operations Plan Documents. EHiring deployment, he stayed in Huntsville; confirmed 
lo^stics and science support arrangements; and maintained communications links with the 
Mission Sdence Teant After deployment, he developed data transmission arrangements with the 
PPs and began compiling and archiving Survey data. 

M. Jarzembski served as a Lidar Scientist for MSFC*s CW CO 2 Doppler lidars. Before 
deployment, he participated in detailed lidar calibration at MSFC and assisted in lidar integration, 
checkout, and cdibration at ARC. During deployment, he participated in routine fidar operations 
and data collection. After deployment, he assisted in lidar demounting and calibration at ARC, 
and in detailed lidar calibration at MSFC. Several months after the GLOBE II mission, he became 
a NASA MSFC employee, where he continued his work with the MSFC lidars. 

D. Cutten served as Aerosol Scientist for two intercomparison flights between the 
Australian Commonwealth Scientific and Industrial Research Organization (CSIRO) F-27 aerosol 
research aircraft and the NASA DC-8 GLOBE aircraft. These flights were conducted near 



Tasmania during GLOBE I. After GLOBE II, he traveled to the United States, joined ESSL, and 
began analyzing the CSIRO and DC-8 intercomparison data. Subsequently, he became System 
Manager and Data Quality Manager for the GLOBE database. 

b. GLOBE Science Coordination D. Bowdle coordinated GLOBE SWG efforts in 
processing and analyzing data from the GLOBE Survey Missions and related GLOBE 
measurement programs. Goals and priorities for this effort were determined fi^om NASA's 
evolving needs and fi-om the sdentific framework developed by ESSL for the GLOBE program. 
Some aspects of this coordination were implemented by routine telephone or mail contact with 
GLOBE investigators, some by annual GLOBE SWG meetings, and some by closer personal 
interaction, as noted below. 

In summer 1989, D. Bowdle and other scientists in the MSEC Aerosol/Lidar group, 
worked closely with Mr. John Porter, Ph D. candidate at the University of Hawaii (UH), during 
his summer internship at MSFC. Mr. Porter analyzed the relationships between the size 
distributions and the visible and infrared optical properties of marine free tropospheric aerosols. 
The data used in this study were obtained by UH researcher Dr. Antony Clarke, during the Fall 
1988 Mauna Loa Aerosol Backscatter Intercomparison Experiment (MABEE) on the Big Island 
of Hawaii. 

In November 1991, D. Bowdle visited the home institutions for most of the GLOBE 
Principal Investigators. The purpose of these site visits was to evaluate GLOBE data processing 
procedures and data quality; to assess selected scientific implications of GLOBE data; to expedite 
release of GLOBE data to the GLOBE database at MSFC; and to provide timely inputs of 
validated GLOBE data to the LAWS development program. 

c. GLOBE Special Sessions D. Bowdle organized special GLOBE sessions in three 
scientific conferences: 

* Fifth Conference, Coherent Laser Radar Technology and Applications, Munich, Germany, June 5-9, 1989; 

* Fall Meeting, American Geophysical Union, December 3-7, 1990, San Francisco, California; 

* Annual Meeting American Meteorological Society, Seventh Symposium, Meteorological Observations and 
Instrumentation, Special Session on Laser Atmospheric Studies, New Orleans, Louisiana, January 14-18, 1991. 

Mr. Bowdle also served as Session Chairman for the Backscatter Session at Munich. 

d- GLOBE Special Issue From September 1987 until March 1991, D. Bowdle served as 
Guest Editor for a special GLOBE section of fifteen papers in the Journal of Geophysical 
Research - Atmospheres (fGR). This effort included soliciting papers, coordinating external 
reviews, performing editorial reviews, and coordinating with JGR. Most editorial responsibilities 
were performed during 1989. The special section was published in March 1991. 

3. GLOBE DATABASE DEVELOPMENT 

a. Hardware Architecture D. Cutten designed and implemented a centralized, remotely 
accessible, interactive workstation architecture for the GLOBE database at NASA MSFC. The 
design centers on a Silicon Graphics 4D25 UNIX workstation (CIRRUS) with magneto-optical 
disk mass storage, ETHERNET connectivity to MSFC GLOBE researchers and other MSFC 
computing resources, and INTERNET connecti\nty to external GLOBE researchers. 
Considerable efibrt was expended in installing, maintaining, and upgrading this vital resource. 

b. Data Archival Software Architecture D. Cutten designed and implemented user- 
fiiendly archival software for the database. Customized software modules read each type of 
GLOBE data, convert it to Hierarchical Data Format (HDF), install standardized file headers, and 



assign standardized filenames. Each module is built on the standard HDF template. HDF, fi'om 
the National Center for Supercomputing Applications, provides several versatile features that 
enhance routine data access for interactive management, retrieval, visualization, and analysis of 
the GLOBE data. In particular, HDF is one of the standard built-in formats for data importation 
into several key data visualization software packages (see next item). Use of HDF also maintains 
compatibility with NASA's Earth Observing System Data and Information System (EOSDIS). 

c. Data Visualization Software Architecture D. Cutten evaluated and implemented 
several powerful object-oriented interactive data visualization software packages for CIRRUS. 
Most data visualization work to date was conducted using the LINKWINDS code from NASA's 
Jet Propulsion Laboratory. Dr. Cutten served as a BETA Tester for LINKWINDS, where he was 
instrumental in prompting numerous strategic improvements in the code. 

d. Data Compilation. Archival, and Distribution Standard data management procedures 
were developed for the GLOBE database. Under these procedures, incoming GLOBE data sets 
were compiled, cataloged, copied onto backup media, and converted to HDF. By the end of the 
contract period, the database had received very little GLOBE I data, but almost all of the GLOBE 
n data except for ancillary meteorological data. Mass data archival of the GLOBE II data was in 
progress. Numerous requests for GLOBE data were received and fulfilled. 

4. GLOBE MODEL DEVELOPMENT 

a. Preprocessor Software D. Cutten, in conjunction with other ES43 Aerosol/Lidar 
Group members, developed several custom pre-processor codes to prepare various types of 
GLOBE data for false-color visualization on CIRRUS. These codes were used to display time vs. 
height cross-sections of aerosol backscatter coefficients from the GLOBE pulsed lidars, and time 
vs. particle size cross-sections of aerosol concentrations from the GLOBE optical particle 
counters. These displays were complemented by line plots of related aerosol and meteorological 
parameters at flight altitude. 

b. Quality Control D. Cutten used the customized data visualization modules to conduct 
detailed studies of GLOBE data quality. These tools revealed several critical, but previously 
undocumented, limitations in the various data sets. These problems were documented and 
corrected where possible, in conjunction with the GLOBE PFs. 

c. Sensor Intercomparison D. Cutten also used the data visualization tools to compare 
contemporaneous data sets fi'om similar GLOBE sensors. These intercomparisons augmented the 
quality control process. Equally importantly, they provided empirical conversion factors fi^om 
commonly measured aerosol parameters to the less commonly measured aerosol backscatter 
coefficients at LAWS design wavelengths. 

d. Aerosol Physicochemical Modeling D. Cutten used measured and analytical aerosol 
size distributions as inputs to Mie scattering codes in computations of aerosol backscatter 
coefficients at 9.1 pm for several common aerosol constituents. The results showed that 
backscatter for mixed particles changes markedly with particle size, type of mixing (internal or 
external), and morphology of mbdng (multiple phases with coated cores, or single phases with 
homogenous mixtures). Dr. Cutten also modeled aerosol backscatter dependence on particle size 
and wavelength for pure ice and water particles. The results of this study were used to interpret 
the empirical conversion factors discussed in item #c above. 

e. Model Synthesis D. Bowdle synthesized results from the GLOBE survey missions, 
other GLOBE measurement programs, and related aerosol data sets, into a model of aerosol 



backscatter properties in the major tropospheric aerosol regimes. This preliminary empirical 
model describes the frequency of occurrence of each aerosol type, and the corresponding 
frequency of occurrence for any given aerosol backscatter value at the 9.1 pm and 2.1 pm laser 
wavelengths. It also describes the atmospheric processes that regulate the regional- and global- 
scale transfer of air mass and aerosol mass from one aerosol domain to another. These model 
studies have confirmed the existence of a remarkably uniform and stable global-scale aerosol 
background population in the middle and upper troposphere. {Roihermel et al., 1989). 

5. AEROSOL LABORATORY DEVELOPMENT 

a. Operation M. Jarzembski, D. Bowdle, and allied MSFC researchers designed and 

implemented an Aerosol Optical Properties Laboratory (AOPL) at MSFC. The AOPL is 
equipped to obtain calibrated measurements of backscatter coefficients from a variety of artificial 
aerosols, using MSFCs 9. 1 1-pm and 10.6-pm CW CO 2 focused Doppler lidars. AOPL research 
objectives include: evaluation of lidar performance theory; development of improved lidar 

backscatter calibration methods; measurement of wavelength-dependent aerosol backscatter; and 
eventual establishment of transfer functions between CO 2 backscatter and aerosol optical 
properties at other wavelengths. M. Jarzembski (an MSFC employee since late 1989) operates 
the AOPL; D. Bowdle provides general scientific oversight. 

b. Measurements M. Jarzembski used the AOPL extensively to calibrate the MSFC CW 
lidars before and after lidar deployment on the GLOBE Survey Missions. Subsequently, he 
conducted exploratory experiments with a new discrete particle cahbration methodology, and 
found excellent agreement with analytical theory for hdar performance in this measurement mode. 
He also compared the performance of the CW lidars using both traditional "hard" target (e.g., 
sandpaper) calibration methodology and the new discrete-particle methodology. 

c. Lidar Measurement Theory D. Bowdle developed a microcomputer-based Monte- 
Carlo simulation code for the CW lidars to explain the details of signal distributions observed 
during the discrete-particle calibration experiments. To complement this code, he also developed 
a simple parameterization for backscatter cross-section distributions in typical populations of 
artificial and natural aerosols. 

6. PUBUCATIONS, PRESENTATIONS, REVIEWS 

a. Publications ESSL researchers published several papers in the open scientific 
literature. D. Bowdle authored two papers and co-authored sbc papery five of these papers were 
included in the GLOBE Special Section of JGR. In addition, M. Jarzwnbski authored one paper 
and co-authored two (one after leaving ESSL), while D. Cutten co-authored one paper. Nine of 
the eleven papers were either directly or peripherally related to GLOBE. These papers are 
itemized in Section UL 

b. Presentations ESSL researchers presented several papers in open scientific 
conferences. D. Bowdle presented nine papers and co-authored four papers in eight conferences. 
S. Wiliams co-authored three papers in three conferences, and D. Cutten co-authored one paper 
in one conference. AD fourteen papers were either directly or peripheraUy related to GLOBE. 
These presentations are itemized in Section in. 

c. Reviews D. Bowdle reviewed one paper for Applied Optics^ two papers for the 
Bulletin of the American Meteorological Society^ three proposals for the proposd authors, and 
three proposals for NASA Headquarters. These re\dews were in addition to the Guest Editor 
reviews for the GLOBE Section of JGR. 
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David Bowdle, June 11-16, 1989: University College, Galway, Ireland. Presented GLOBE results 
at the International Conference on Aerosols and Background Pollution. 

Steve Williams, June 10-14, 1989; NASA Ames Research Center, Mountain View, California. 
Reviewed plans for the First GLOBE Survey Misssion. 

Steve Williams, July 11-29, 1989: Auckland, Wellington, and Christchurch, New Zealand; 
Melbourne, Canberra, and Sydney, Australia; Tokyo, Japan; Honolulu, Hawaii; and San 
Francisco, California. Conducted site evaluations for the First GLOBE Survey Mission. 

Steve Williams, August 28-29, 1989: NASA Headquarters, Washington, D.C. Reviewed plans for 
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Meering Activity Planning and Support 
Prepared by: Laurie Collins 


Planning, logistics, administrative support during the meeting and proceedings 
preparation were undertaken for the following meetings. Facilitation of attendance 
for the Science Team was provided for workshop participants as needed and travel 
support was provided as necessary to the Science Team when attending other 
meetings as a designated member of the lAWS Science Team. 


Meeting Activity: 

February 1989 - New Orleans, lA - Geostationaiy Platform Science Steering Committee 

March 1989 - MSFC, AL - ESGP Phase A Meeting 

May 1989 - Huntsville, AL - GLOBE Science Working Group Meeting 

May 1989 - Washington, DC - Geostationary Platform Meeting 

July 1989 - Huntsville, AL - LAWS Science Team Meeting 

October 1989 - MSFC, AL - ESGP Phase A 

January 1990 - Pasadena, CA - LAWS Science Team Meeting 

August 1990 - Boulder, CO - lAWS Science Team Meeting 

February 1991 - Clearwater, FL- lAWS Science Team Meeting 

July 1991 - Snowmass- Aspen, CO - LAWS Science Team Meeting 

December 1991 - Toulousse, FRANCE - lAWS Science Team Meeting 

December 1991 - Huntsville, AL - lAWS Science Team Meeting 

December 1991 - Guemeville, CA - US Science Working Group Meeting 

July 1992 - Falmouth, MA - LAWS Science Team Meeting 



A Regional Scale Pre-Convective Simulation Incorporating A Satellite- 

Derived Data Assimilation 

By 

J. Aaron Song, Richard T. McNider and Dan Casey 


1. Introduction 

Supported by NASA Contract No. NAS8-37585 we have been making initial computations to 
determine whether high resolution boundary layer modeling in conjunction with satellite derived moisture 
fields can explain observed meso-bcta convective patterns. Mesoscale numerical modeling experiments have 
been performed on a selected pre-convective environment in which synoptic -scale dynamic forcing is small, 
while localized convective developments are observed under a favorable condition. The favorable condition 
seems to be best realized when upward velocities in a vertical motion field generated by the local terrain 
effect overlaps with higher moisture in a satellite-observed lower-tropospheric moisture field. In this report, 
we first summarize the numerical experiments that have been performed in which diagnostic analysis shows 
observed storm developed over the locations which arc predictable using the GOES satellite data and a 
mesoscale boundary layer model with high resolution resolving the complex terrain effects. Then, we 
introduce a newly-derived technique which is designed to prognostically incorporate the satellite hourly 
surface observation into a mesoscale boundary layer model such that the satellite information is assimilated 
in a dynamically consistent manner. 

2. Diagnostic Mesoscale-Satellite Coupled Analysis 
2.1 A Moist-Thermal Mechanism 

One of the essential mechanisms by which cumulus clouds develop over a synoptically-undisturbed 
quiescent environment with complex terrain is the breaking of the coping inversion through a positive 
feedback effect involving a local thickening of the mixed layer and the itKloud latent heat release. Figure 1 
shows a schematic illustration of a possible breakthrough of the capping inversion by cumulus 
development within the boundary layer to which terrain-induced mesoscale lifting contributes. If the local 
lifting condensation level (LCL) is above the mixed layer, no clouds form, but if the LCL is within the 
mixed layer, dry thermals in the boundary layer, with turbulent velocities generally on the wder of 1ms 
can produce shallow cumuli within the layer between the LCL and the inversion base. The cumuli may not 
be strong enough to penetrate upward if there is no supportive forcing on a larger scale. In the case when 
terrain-induced mesoscale lifting exists, as shown in figure 1. the mesoscale upward motion produces both 
a moisture convergence (which lowers the LCL) and a rise of the mixed layer top (therefore thickens the 
‘ layer of the moist thermals). The thickness of the moist thermals is dynamically important, since, once 
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thickened, the parcels within the cumuli would have a vertically longer path to strengthen, resulting in a 
larger latent heating and a possibly stronger overshooting at the cloud top. This results in a thicker mixed 
layer and, consequently, stronger turbulent velocities. The stronger turbulence then, as supported by the 
continuous mesoscale lifting, produces stronger dry thermals below the LCL and moist thermals above the 
LCL, and the above process is repeated. Once this positive feedback is established, it can result in a 
breakthrough of the inversion and deep cumulus convection may initiate thereafter. 

2.2 Mesoscale-Satellite Coupling 

Mesoscale modeling experiments have been performed over the selected domain shown in figure 
2a, in which the complex terrain surrounding the Smoky Mountains is contoured on a lOkm-by-lOkm grid. 
This domain is chosen from considerations of (1) the GOES VISSR visible and IR (with pixel resolution 
upto Ikmxlkm), as well as the VAS, data must in general be averaged to be used, avoiding the error 
associated with the solar zenith angle (Gautier et al.,1980); (2) covering enough mountainous area in order 
to investigate the generation of vertical motions by complex terrain; and (3) use of in-situ observation for 
model verification from the data of COHMEX ( Williams et al.,1987). The mesoscale model produced 
vertical motion field on the lOkm-by-lOkm grid (the simulation is initialed with the 1400UTC sounding of 
19 June,1986, of COHMEX) is shown in figure 2b. Combinning figures 2a and 2b, it can be understood 
by the good corelation between the terrain and the vertical motions that it is essnetial to use a high- 
lesolution boundary layer model which is able to at least resolve the mesoscale features of the complex 
terrain in this cae. A 20km-by-20km grid has also been tested (not shown) which does not result in the 
vertical motions essential to the study purpose of this research. 

A split-window retrieval technique (Jedlovec, 1990) was used on the 19 June, 1986, case using the 
multi-spectral VAS (the window<hannel and the lower-tropospheric moisture channel) data which results in 
a mesoscale distribution of the vertically-integrated precipitabic water around 1600UTC over the COHMEX 
domain. Because different domains are used between this research and that of the retrieval, a Barnes 
objective analysis is used to analyze the precipitable water data onto the domain of this study, shown in 
figure 3. As limited by the retrieval data coverage, the precipitable water data ovct the eastern third of the 
domain shown in figure 3 (which, unfortunately, includes the Smoky Mountains) must be used with greater 
care as compared with the rest of the domain. Nevertheless, a moist-tongue can be identified ftom the 
precipitable water data which extends southwestwardly from the southeast Kentucky, across Tennessee to 
northern Alabama, with a local moisture peak around the central Tennessee. 

A diagnostic analysis was made which illustrates the mesoscale distribution of the thickness of 
moist thermals as produced using the model PBL top height and the LCL obtained from a moisture 
enhencement scheme. The moisUue enhencement scheme is designed to incorporate the satellite moisture 
pattern into mesoscale modeling in which the moisture field is otherwise a completely passive model 
variable. Currently, the scheme uses a base surface moisture field (a constant value), and multiply to it 
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grid-by-grid a ratio defined as dividing the grid-point prccipitable water by the domain-averaged precipitaWc 
water. This results in a more realistic surface moisture distribution over the mesoscale domain, which 
would not be possible to be generated from using simply an ordinary mesoscale boundary model such as 
that used in this research (originally developed by Pielkc.1974). With this refined surface moisture, the 
LCL can be calculated, and, by subtracting the LCL from the model PBL height, the moist thermal 
thickness field is obtained, shown in figure 4. From figure 4, it is seen that over the central Tennessee and 
over the Smoky Mountains, larger positive values of the moist thermal thickness arc found, whereas 
negative values can be seen surrounding the Smoky Mountains, Figure 5 shows the GOES VISSR visible 
image of 2(XX)UTC, 19 June, 1986, on the lOkm-by-lOkm grid, whch indicates storm developments, as 
seen from the much larger brightness values over the cloud tops, in cerlral Tennessee and over the Smoky 
Mountains. Qualitatively, this meso-beia scale convective pattern seems to be predictable, or at least 
explainable, using the aforementioned mesoscale boundary layer model in conjunction with the satellite- 
derived moisture enhancement scheme. 

3. On-Going Research 

Incorporation of satellite information in regional scale atmospheric modeling may be denoted as a 
numerical modeling with one type of Four-Dimensional Data Assimilation (FDDA). Generally, on the 
global scale the FDDA has been carried out during the past two decades with the use of the balanced- 
equation constraint through the geostrophic adjustment process. On the mesoscales, or scales smaller than 
the Rossby radius of deformation, however, the FDDA seems currently subject to a certain degree of 
uncertainty. The problem seems to not just concern that the balanced-equation may not hold on 
mesoscales, but also that the traditionally linearly-derived radius of deformation may not be itself 
meaningful (Stauffer et al.,1991). A single-level (surface) data assimilation does not seem to result in 
improvement on the regional scale simulations (Dr. N. Seaman; personal communication). Even with a 
deep-layer information, for instance using the VAS data, Hayden and Schmil (1991) indicated that without a 
proper dynamic constraint the current satellite-derived FDDA does not seem to be able to consistently 
improve on the mesoscale numerical modeling. 

An experimental thermodynamic/dynamic constraint linking the GOES hourly data and mesoscale 
boundary layer modeling is currently under development in this research. Figure 6 shows schematically the 
computational steps of assimilating satellite skin temperature information (in terms of hourly changing 
rates) into the mesoscale PBL modeling which determines the surface temperature and moisture as well as 
the sensible and latent heat fluxes across the air-ground interface. In McNider et al.(1993) this data 
assimilation approach is introduced and discussed using the one-dimensional data from FIFE (First ISLSCP 
Field Experiment; where ISLSCP: IntOTiational Satellite Land Surface Climatology Program), while in 
Song et al.(1993) three-dimensional experiments incorporating this data assimilation will be discussed with 
real-case simulations using observations from COHMEX and CaPE ( Williams el al.,1992). 
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HGURE CAPTIONS 


Figure 1: A schematic illustrauon showing a more effecUve breakthrough of a capping inversion by cumuli 
development between the LCL and the PBL top under a mesoscale terrain-induced lifUng, in contrast to a 
less-likely breakthrough in which no mesoscale lifting contributes. 

Figure 2: (a) Topography heights (m) plotted on a 4.5km-by-4.5km grid over the study domain, with 
contouring interval of 200m; and (b) Model produced vertical motions (cm/s), averaged over the layer of 1- 
2km above the surface, on a lOkm-by-lOkm grid with solid/dashed lines indicaUng upward/downward 
velocities and a contouring interval of 10 cm/s. 


Figure 3: Precipitable water (cm) distribution on a lOkm-by-lOkm grid over the study domain at 1600UTC, 
19 June, 1986, as obtained from the use of a split-window retrieval technique (developed by Jedlovec,1990). 

Figure 4: Mesoscale distribution of the thickness (m) of Moist Thermals (defined as the layer thickness 
between LCL and PBL top) around 1800UTC, 19 June.1986. on a lOkm-by-lOkm grid, obtained from a 
mesoscale-satellite coupled diagnostic analysis (see text for details). 

Figure 5: The GOES VISSR visible image in brightness (0-255). as plotted on a lOkm-by-lOkm grid, on 
2000UTC, 19 June.1986, over the study domain. Relatively much higher brightness (around 2(X)) show 
the cloud tops of deep convection, while lower brightness (80-100) generally indicate surface under dear- 

sky. 

Figure 6: A schematic illustration showing the computational steps of an experimental data assimilation 
technique developed in this research, which shows the incorporation of satellite skin temperature hourly 
changing rates into a prognostic surface energy balance equation in such a way that the assimilated satellite 
infemnation is dynamically and thermodynamically consistent with that of the mesoscale model (details are 
discussed in McNider et al.,1993 and in Song et al.,1993). 
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Prepared by 


Douglas Mach 


I. Task 1 

High Altitude Aircraft Instrument Package 

My part of task 1 was to help with the integration of the presently available 
instruments onto the new platform. I also participated in the selection of new in- 
struments and experiments to be built and integrated into the new platform. 

II. Task 2 

Optical Pulse Analysis 

This task dealt with the analysis of presently available optical pulse data col- 
lected on the many U-2 storm overflights made by the lightning group at 
NASA/MSFC. The data set consisted of a huge number of pulses taken from many 
storms over several years. Because of the immense number of pulses in the data set, 
it was ideal for performing statistical pulse analysis. My major effort in this sub-task 
was concentrated on perfecting an automatic program to find and categorize the op- 
tical pulse data. I was able to correct and debug a program to do automatic analysis 
of the pulse data and I have tested the program on smaller subsets of the data base. 
My task in the next year was also to have the program do the analysis on the com- 
plete data set and to determine the relevant pulse characteristics. These included 
pulse rise times, pulse duration, pulse shape, pulse fall time, and optical and electric 
field arrival time delays. Once the statistics were determined, I compared my results 
to those of others to determine any differences. 

Cloud Light Scattering Model 

One of the questions raised in previous work in the detection of optical 
pulses from clouds is that the optical signals seem to lag the electric field pulses. In 
addition, the optical pulses detected from above the cloud have different charac- 
teristics than the pulses from lightning detected below the cloud. One possible ex- 
planation for this is that the scattering and absorption of light by cloud particles 
delays and distorts the optical pulse as it passes through the cloud. I tested this 
theory by modeling the interaction of the cloud particles and the photons by using a 
simple Monte Carlo cloud model to simulate the interaction of cloud particles with 
photons. This program yielded some interesting results. I produced a large data 
base of pulse delay, and angular, and spectral dependence for the spectral bands 
that are detected by the NASA/MSFC optical instruments. The results of the cloud 
model studies were then used as calibrations for theoretical studies of the optical 
scattering process. 



